Measurements of Accelerations of Large Neutrally-buoyant Particles in Intense 

Turbulence 
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We measure acceleration statistics of neutrally buoyant spherical particles with diameter 0.4 < 
d/i] < 27 in intense turbulence (400 < Rx < 815). High speed cameras image polystyrene tracer 
particles in a flow between counter- rotating disks. The measurements of acceleration variance, (a^), 
clearly resolve the transition from the tracer like behavior of small particles to the much smaller 
accelerations of large particles. For d > 5ri, (a^) decreases with diameter as d~^^^ in agreement with 
inertial range scaling arguments. A model relating (a^) to the pressure structure functions matches 
the transition from small to large particle behavior if the particles respond to pressure differences 
over (1.7±0.3)ci. A model relating {a'^) to the fluid acceleration averaged over the particle diameter 
predicts the transition with no free parameters, but does not show clean inertial range scaling in 
the size range studied. Consistent with earlier work, we find that the scaled acceleration probability 
density function shows very little dependence on particle size. 
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The study of accelerations of turbulent fluids and of 
particles entrained in turbulent fluids has produced many 
recent insights [ij. Even excluding differences between 
turbulent flows, the parameter space to describe parti- 
cle motion contains a wide range of possibilities span- 
ning both particle size and particle density. The most 
intensely studied region of this parameter space has been 
small particles that are heavier than the fluid 0, 13| • The 
motion of small heavy particles is relevant to many nat- 
ural systems including water droplets in clouds and sed- 
imentation. This problem is also attractive because it is 
accessible numerically by modeling the particles aspoints 
interacting with the fluid through a drag law 0, Sig- 
nificant study has also focused on the region of the parti- 
cle parameter space in which the particles have positive 
buoyancy, such as air bubbles in water 0]. Our in- 
terest in this paper is in the accelerations of neutrally 
buoyant particles that vary in diameter from tracer par- 
ticles {d = 0.477) to much larger than the Kolmogorov 
length (d — 27 rf). A number of applications, including 
many marine organisms, inhabit the large neutrally buoy- 
ant domain. Additionally, tracer particles used in liquid 
phase fluid measurements are often neutrally buoyant. 

Several previous experiments have addressed the prob- 
lem of large neutrally buoyant particles in turbulent flow. 
Voth et aZ.(2002) measured size dependence of particle 
accelerations in order to confirm tracer particle behav- 
ior in their experiment. They have relatively few data 
points; however, they suggested that the accelerations of 
these large particles might be modeled by terminating the 
cascade at the particle diameter, d. The particle acceler- 
ation then scales with the fluid acceleration at the scale 



of the particle diameter, which implies (a 
in the Kolmogorov inertial range. Qureshi et 
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sured accelerations of neutrally buoyant helium bubbles 
in a wind tunnel. Their particles ranged from 7?] to 
2577 = Q.IL which at their relatively low Reynolds num- 
ber {R\ = 160), spans from the dissipation scales to the 
injection scales. But they are unable to observe the tran- 
sition from dissipation to inertial range. They highlight 
the fact that the particle acceleration variance can be 
connected to the pressure structure function evaluated 
at the diameter of the particle, (aparticio) °^ (AP^(d))/(i^ 
which gives the same d"^/"^ scaling in the inertial range 
for pressure structure functions that have r'*/"^ inertial 
range scaling. 

Numerical simulations of this problem are difficult 
because a complete model must resolve the turbulent 
scales and the boundary layer around individual par- 
ticles [1, [1^. Calzavarini et al. overcome the severe 
Reynolds number limitations of full simulations by using 
effective equations for particle motion that utilize phe- 
nomenological Faxen corrections [llj . 

In this letter we present experimental measurements 
of the acceleration of neutrally buoyant particles at large 
Reynolds numbers and a range of particle size that 
spans from dissipative lengthscales into the inertial range. 
These measurements are the first to clearly resolve the 
transition from tracer particle to large particle behavior 
and allow a stringent test of various models that have 
been proposed to quantify the accelerations of large par- 
ticles. 



Lagrangian particle tracking measurements were carried 
out in the von Karman fiow between counter rotating 
disks described in The relevant parameters of the 

flow, including Rx and the Kolmogorov microscales, are 
shown in Table |T1 Three Phantom v7 cameras operat- 
ing at 20,000 frames per second and 256 x 256 pixels per 
frame were arranged in the central horizontal plane and 
focused on a 1 cm'^ region at the center of the flow. A 
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TABLE I: Table of flow parameters: /, frequency of the rotating disks; u, kinematic viscosity of the fluid; Rx, Taylor Reynolds 
number; u, rms velocity of the flow; L, energy input length scale; rj, Kolmogorov length scale; r,,, Kolmogorov time scale; Nf, 
number of frames in each Kolmogorov time; imaging volume; and Aa;, the distance in the flow corresponding to one pixel. 



frequency-doubled, pulsed Nd:YAG laser with 50W av- 
erage power was used for illumination. Two beams were 
required to allow for forward scattering in all three cam- 
eras. Vertical polarization minimizes secondary reflec- 
tions from the particles [13]. 

All particles used were polystyrene spheres {p — 1.05 
g/cm'^). Measurements were made with the flow seeded 
with both mono-dispersed particles from Duke Scientific 
with known diameters (d=26, 55, 134, 222, 300, 400 ^m) 
and with poly-dispersed 'grinding media' with diameters 
in the range 600-990 /im obtained from Norstone, Inc. 
We use a reflection pair method to determine the diam- 
eter of the large particles from the images of each tra- 
jectory. The central camera records two reflections from 
each particle in the imaging volume because of the two 
beams. For large particles {d > 100 ^m) the two reflec- 
tions can be resolved, allowing measurements of particle 
size in every frame. Excluding particles whose diame- 
ter changes over time minimizes error due to the small 
fraction of non-spherical particles in the grinding media. 

An 8% NaCl solution by mass {v = 1.29 x lO^m^/s, p = 
1.048g/cm'^) is used in order to density match particles 
with diameters greater than 100 microns, but deionized 
water is used for particles with diameters of 26 and 55 
microns. The effect of a 5% density mismatch on par- 
ticle accelerations is nearly negligible for all our particle 
sizes however keeping the large particles suspended 
in the fluid requires careful density matching. 

We use the data acquisition system used in Im- 
ages are acquired in one second sequences, which take 
3-5 minutes to download. Particle positions extracted 
from the images are stored and later tracked. For each 
particle size, we obtain on average 65 million particle 
positions from 40-500 sequences. Accelerations are de- 
termined with a quadratic fit to the particle trajectories. 
Figure [T] shows the measured acceleration variance for 
one data set as a function of fit time. For the highly 
intermittent accelerations of particles in turbulence, it is 
very difficult to obtain fits that are not affected either by 
measurement errors or by smoothing of the trajectories. 
We use an established method [l2| to extrapolate back 
to zero fit time in order to determine the acceleration 
variance. 

All results reported are for analysis of two dimensional 
trajectories. With three cameras, it is possible to deter- 
mine 3D trajectories. However, we encountered difficul- 



ties with stereomatching from images of particles that 
fill a large fraction of the field of view and so the 3D 
tracks are extremely fractured. As a result, we find much 
better statistical convergence for the 2D trajectories. In 
this letter, we focus on the axial component of accelera- 
tion (parallel to the axis of rotation of the disks in our 
flow). The transverse component shows similar results 
with larger measurement uncertainties due to the illumi- 
nation geometry. In the Reynolds number range studied, 
the transverse acceleration variance of fluid particles in 
this flow is 10% to 20% larger than the axial variance p^ . 

Figure [2I^a) shows acceleration variance measurements 
as a function of particle diameter. The acceleration is 
normalized by Kolmogorov variables to give qq. For small 
particle diameters {d/rj < 5) the normalized accelera- 
tion variance is nearly constant. At larger diameters 
the normalized acceleration variance falls off implying 
that these large particles do not simply trace the flow. 
There is no detectable dependence on Reynolds number. 
However, at the lower Reynolds number, the largest par- 
ticles studied are only 10.97?. For fluid particles, qq is 
found to vary by less than 17% over the Reynolds num- 
ber range we study [T^ . Data is presented for two meth- 
ods of determining particle positions. The first method 
treats each reflection observed in an image as a parti- 
cle. This method works well for mono-dispersed particles, 
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FIG. 1: An example extrapolation of acceleration variance 
back to zero fit time. Acceleration variance measurements for 
single reflection measurements of the 300/im {d/rj = 10.64) 
particles taken at Rx = 717 are plotted at various fit times, 
along with the statistical error in these values. 
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d < 400/Ltm, and we can use data from all three cameras 
to improve statistical convergence. The second method 
was developed to determine the size of polydisperse par- 
ticles from the images. Here we use only the central 
camera, and require two reflections to be observed that 
move parallel to each other. This analysis can be done 
for particles with diameter d > 134/xm. Data points are 
plotted for both of these methods to show that no sys- 
tematic error is introduced when we use the reflection 
pair method. Error bars in Figure [2] represent statistical 
uncertainties. Systematic errors associated with determi- 
nation of the energy dissipation rate and extrapolation 
also exist [l^l • Figure [2l^b) includes data from two pre- 
vious measurements of accelerations of large particles. 

Figure ^c) shows the measured acceleration variance 
compared with predictions from three models. The first 
of these models is the d~^^^ scaling predicted for particles 
with diameters in the inertial range . Our data shows 
good agreement with the d^^^'^ scaling for d > brj. 

A second model is suggested by the connection be- 
tween the particle acceleration variance and the pressure 
structure function 0]. In addition to showing d"^/^ iner- 
tial range scaling, this Pressure Difference model can be 
extended to predict the acceleration variance for small 
particle diameters d < ij using the pressure structure 
function for separation distances in the dissipation range. 
In a previous experiment in the same apparatus that we 
use, Xu et al. [15| measured the pressure structure func- 
tion over the range O.Tr; < r < 140?7. We use their data 
to create the prediction shown as a dashed line in Fig. [21 
When the separation at which the pressure difference is 
sampled is equal to the particle diameter, this model pre- 
dicts larger accelerations than we measure at large parti- 
cle sizes. However, this model is a good fit when we allow 
the separation to be an adjustable parameter and sam- 
ple the pressure structure functions at (1.7±0.3)d. This 
raises questions about the mechanism which would cause 
the particle's acceleration to be determined by length- 
scales larger than its diameter. 

The third model assumes that the particle acceleration 
is equal to the spatial average of the fluid acceleration 
over the volume of a particle. This model closely matches 
the model used by Calzavarini et al [llj. Their model 
differs by a drag term that contributes only 1% to the 
neutrally buoyant particle acceleration. 

In this Fluid Acceleration Averaging model, the parti- 
cle acceleration variance is given by 

(aparticic') = / / {af{v)af{v'))d^vd^v'. (1) 

J V ^ v' 

Integrating and using the isotropic expression, Rij = 
RNN^ij + {Rll — RNN)rirj/r'^ for the acceleration cor- 
relation function yields 
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Xu et al. [l5| have measured the acceleration correla- 
tion functions, RLhij") and RNN{f), in the flow that we 



FIG. 2: (color online) (a) Acceleration variance as a function 
of particle diameter. Unfilled symbols represent data taken 
at the lower disk frequency and hence lower Reynolds num- 
bers {R\ — 396,449). Filled symbols represent data taken at 
the higher disk frequency and hence higher Reynolds num- 
bers {Rx = 717,813). Triangles are single reflection mea- 
surements. Diamonds are reflection pair measurements, (b) 
Filled and unflUed symbols are the same as in A, but pair mea- 
surements are not shown for data sets where single reflection 
measurements are available. X's are measurements from [l^ . 
Grey squares from |j8|]. The solid curve is the the inertial scal- 
ing derived from the cascade truncation model, (c) Symbols 
are the same as in A. The solid curve is the d~^^^ inertial 
range scaling. The dashed curves show the predictions of the 
Pressure Difference model. The red dashed line uses a sepa- 
ration equal to the particle diameter. The dashed orange line 
uses the separation which best fits fit to our data (1.7d). The 
dotted black curve is the prediction of the Fluid Acceleration 
Averaging model. 



study. Integrating their measured functions in Eq. [5] pro- 
duces the dotted curve shown in Fig. [2][c) . This model 
gives a very good fit to our data both at small scales 
and in the transition to inertial range scaling. No free 
parameters are needed to adjust the length scale in the 
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a/<a2>i/2 

FIG. 3: Probability density functions of acceleration at Tay- 
lor Reynolds numbers (a) 717 and (b) 396. Particle sizes in 
d/f) represented by each symbol are as follows: black unfilled 
diamonds, (a) 1.95 and (b) 4.75; blue filled diamonds, (a) 3.23 
and (b) 7.87; green unfilled triangles, (a) 4.37 and (b) 10.64; 
and red filled triangles, (a) 5.82 and (b) 14.18. 



model to the measurements as is required for the pres- 
sure structure function model. However, this third model 
predicts accelerations that are somewhat too small for 
large particles sizes. It also does not display c?"^/"^ scal- 
ing in the range of particle sizes we studied. The accel- 
eration correlation functions do not show inertial range 
scaling until d > iOrj 15], and so inertial range scal- 
ing for this model is not expected in the range d < 25?7 
that we study. Although the data seems to match the 
d~'^^^ scaling quite well, the lack of scaling in the accel- 
eration correlation functions suggests the possibility that 
the observed d~'^/^ scaling for small particles sizes is only 



approximate and that much larger particles sizes (along 
with high Reynolds numbers) are needed before rigorous 
d^"^/^ scaling will exist. Figure [3] shows the probability 
density function (pdf) of particle accelerations normal- 
ized by the standard deviation for particle sizes d < I5r]. 
Within our measurement error, these pdfs agree with the 
form of the fluid acceleration pdf In the tails of 

the pdfs, the probability densities for large particles sizes 
may be slightly below the fluid particle pdf. This is more 
pronounced in the data from the large poly-dispersed par- 
ticles (not shown). However, in looking carefully at our 
data we found errors due to particle finding, tracking, 
and nonspherical or stray particles that seem to affect 
the rare events by as much as any discrepancy with the 
fluid acceleration pdf, so we do not draw any conclu- 
sions from these slight differences. This particle size in- 
dependence of the acceleration pdfs has been observed 
before 0, [l^ . flTj and remains somewhat puzzling. 

These measurements of accelerations of spherical neu- 
trally buoyant particles in intense turbulence clearly re- 
solve the transition from the tracer regime to the large 
particle regime. Our measurements of the acceleration 
variance as a function of particle size show the best agree- 
ment yet measured with the inertial range scaling. 
We present two models that capture the transition from 
tracer to large particle behavior. The success of these 
simple models is perhaps surprising given the complex 
flow in the boundary layer around large particles in tur- 
bulence. The parameter space of this problem is vast, 
but these measurements seem to put our understanding 
of one region, large neutrally buoyant particles, on a solid 
foundation. 
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